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In order to switch its direction, a macroscopic magnetic moment of a superparamagnetic system 
has to overcome an energy barrier which is due to magnetic anisotropies in the material. The 
switching rate is usually described by a Neel-Brown-Arrhenius law, in which the energy barrier 
is a constant for a given temperature. However, magnetic anisotropies are in general depending 
on temperature. We take this temperature dependence into account and derive a generalized Neel- 
Brown theory for the switching rate. Importantly enough, this generalization has a tremendous effect 
on the physical interpretation of switching processes obtained in experiments on the basis of the 
conventional Arrhenius law. In particular, the assumption of a temperature-independent anisotropy 
may result in an overestimation of the attempt frequency by several orders of magnitude. 


PACS numbers: 

The development of devices for magnetic storage media 
faces several challenges in the ongoing miniaturization 
of information units. One of the fundamental physical 
problems is the so-called superparamagnetic limit. In the 
small scale limit the alignment of the ferromagnetic par¬ 
ticle with a certain axis is no longer stable but affected by 
thermal energy. Despite similarities between para- and 
superparamagnetism, distinct differences exist. Quan¬ 
tum mechanics allows to describe paramagnetism rigor¬ 
ously with the magnetic field as a small perturbation. 
In superparamagnetism, the spins of a nanostructure are 
exchange coupled which causes a quasi-classical behavior 
with a large variety of possible energetically continuous 
states. Due to the large number of coupled spins, infe¬ 
rior energy contributions of the single spins become im¬ 
portant (dipolar, spin-orbit effects). Such effects summa¬ 
rized as anisotropies compete with the thermal energy in 
a system of small spatial dimensions. The total energy 
of the nanomagnet reveals remarkable differences com¬ 
pared to the single moment case. In particular due to the 
competition with thermal energies, strongly changing be¬ 
havior appears as function of temperature, for instance, 
blocking of the macrospin along a certain direction (easy 
axis of magnetization). In the latter situation, an en¬ 
ergy barrier between the two degenerate states oriented 
parallel to the easy axis suppresses thermal switching. 
In case of an uniaxial system, the barrier height scales 
with the volume and strength of the anisotropy. At large 
enough temperature, the barrier can be overcome and 
the moment can flip back and forth. The motion of the 
magnetic moment in the potential well allows to define 
the so-called attempt frequency of the reversal process 
[1-3]. Besides a Boltzmann factor, scaling with the ra¬ 
tio of the total anisotropy energy and the thermal en¬ 
ergy, the attempt frequency is in general considered as 
the only quantity which is influenced by temperature. In 
general, experiments are explained under the prerequi¬ 
site of constant barrier heights to determine the mag¬ 
netic anisotropy, thereby largely ignoring the extremely 
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FIG. I. (Color online) Sketch of magnetic potentials of su¬ 
perparamagnetic nanoislands for different temperatures with 
easy axis in the out-of-plane direction. To switch between the 
minima, the magnetization M has to overcome the anisotropy 
energy barrier which typically depends on M(T). With in¬ 
creasing temperature, both the magnetization and the energy 
barrier, are decreased. 

varying attempt frequencies that are simultaneously ob¬ 
tained [4]. That unphysical results may arise has been 
pointed out in Ref. [5]. The authors have realized that 
the anisotropy may be temperature dependent and thus 
can lead to a wrong value of the attempt frequency. In 
general, it is well known that the anisotropy changes with 
temperature due to changes of microscopic lattice param¬ 
eters, strain and others. In this Letter, we point out that 
besides the blocking temperature T^, a second tempera¬ 
ture scale becomes relevant in the macrospin description 
which is the ordering temperature Tc- So far, this fact 
has not be taken into account for the interpretation of 
experiments with superparamagnets. We show that for 
an appropriate interpretation of the temperature behav¬ 
ior of nanoparticles, temperature has to be considered 
with respect to the ordering temperature. In particu¬ 
lar, we demonstrate that the finite ordering temperature 
has a strong impact on the switching dynamics of super¬ 
paramagnets. We illustrate that experimentally deter¬ 
mined prefactors [6-9] have to be reinterpreted particu¬ 
larly when the temperature scale of the measurements is 
comparable with the ordering temperature. 

Microscopically a superparamagnetic system can be de- 
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scribed by the Heisenberg model 

= -J^Si-S,+i/aniso, (1) 

(«d> 

where the spins tend to align parallel due to the ex¬ 
change interaction of strength J in order to form a macro¬ 
scopic magnetic moment. The exchange interaction itself 
is rotationally invariant and thus does not favor a certain 
direction for the resulting magnetic moment. However, 
superparamagnetic systems typically possess either in¬ 
trinsic material symmetries or external ones such as a 
shape anisotropy. Then, distinct directions of the mag¬ 
netic moment are preferred. To model this feature, a 
generic anisotropy term is introduced in Eq. (1). The 
specific form of the anisotropy has to be defined for the 
particular experimental situation. 

We assume the magnetic system as large enough to 
form a macroscopic spin but still sufficiently small to 
avoid a separation into multiple magnetic domains. In 
this superparamagnetic limit, we can replace the collec¬ 
tion of microscopic quantum spins by classical vectors 
and assume a translationally invariant system 
((S'?) = (S^)Vi). Changing to magnetization as order 
parameter the uniaxial anisotropy is described by an en¬ 
ergy density and the anisotropy part in Eq. (1) becomes 

-ffaniso = -KV COS^ 6, (2) 



FIG. 2. (Color online) Temperature dependence of the 
magnetization m — M(T)/M(0) (full lines) and its square 
m? (dashed lines) described by (a) the Brillouin function 
S(S, Tc) for a spin S = (b) Langevin function L{Tc) for a 

classical spin, and (c) a linear dependence m{T) = 1 — TjTc- 
Tc denotes the Curie temperature of the system. 


with AE(T) typically decreasing with increasing temper¬ 
ature, as it is sketched in Eig. 1. 

We focus in this work on the consequences of a 
temperature-dependent anisotropy for the switching rate 
and thus investigate a model where the rate is given by 
a generalized Neel-Brown-Arrhenius law 


KnM(T)2 

f{T) = foe 


( 4 ) 


with 0 being the angle which the magnetization encloses 
with the z-axis. The total Hamiltonian for a parti¬ 
cle of volume V and M has two energy minima for 
M = ±|M|e; 2 . In order to reverse its direction, the 
total magnetic moment has to switch from one mini¬ 
mum to the other by surmounting the energy barrier 
AE = H,niso{0 = 90°) - ffaniso(0 = 0°, 180°) = KV 
while crossing the highest energy state at ^ = 90°. The 
interaction part of Eq. (1) is irrelevant for determining 
the energy barrier as it is rotationally invariant and will 
only be considered implicitly as the mechanism causing 
the magnetic order. At finite temperature, the difference 
of the free energy AF{T) = AE{T) - TAS{T) deter¬ 
mines the height of the switching barrier. However, we 
assume the entropy S{T) does not change significantly 
by a rotation of the magnetization and set TAS <C AE. 
Eor the discussion in the following, it is crucial to note 
that the anisotropy is related to the particle magnetiza¬ 
tion via a conventional power law K = K[M(T)/M{0)]^ 
with the temperature dependent magnetization M(T) 
and M(0), K the corresponding values at zero Kelvin 
[10-14]. Eor the sake of simplicity we consider the uni¬ 
axial anisotropy to be of dipolar origin. Then, the en¬ 
ergy difference between hard and easy axis is given by 
AE = VNM‘^[lb] with N stemming from the demagne¬ 
tization tensor. Thus, we find 

AE{T) = VNM(Tf = KoM{Tf , (3) 


We fix the attempt frequency to /o = 10^ s“^ and the 
anisotropy energy to KoM{T = 0)^ = 1 eV throughout 
this work. The conclusions below hold even for reason¬ 
able deviations from Eq. (4). Eor example, a different 
exponent u in M{T)^ may be absorbed in the actual tem¬ 
perature dependence of M(T)^ while a (non-exponential) 
temperature dependence of the prefactor /o [2] is in first 
order negligible. 

The actual magnetization of a physical system is af¬ 
fected by a multitude of material properties, such as in¬ 
trinsic interactions or the finite sample size. While in 
ferromagnetic bulk samples the temperature dependence 
of the magnetization is typically well described by a Bril¬ 
louin function, this is often not true in systems of reduced 
dimensions [13, 16, 17] . Thus, a general assumption 
of M(T) cannot be made without explicit knowledge of 
the system. To illustrate that the explicit temperature- 
dependence of M(T) has a major effect on the switching 
rate, we consider three different examples in this Letter: 
(a) the Brillouin function describes quantum spins in a 
bulk sample, (b) the Langevin function refers to a clas¬ 
sical moment (M ^ oo), and (c) a linear dependence 
such that M(T) = M(0)(1 — TjTc) (with the Curie tem¬ 
perature Tc) as a possible variant of an unconventional 
magnetization curve which can be realized in magnetic 
nano-islands [11]. These functions and their squares are 
shown in Eig. 2. 

In general, M(T) decreases with increasing temperature 
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FIG. 3. (Color online) Decadic logarithm of the switching 
rate Eq. (4) for the temperature-dependent magnetizations 
shown in Fig. 2. Strong deviations from an Arrhenius law 
with a constant anisotropy KqM{TY = i^oM(O)^ (dashed 
line) appear even for T <C Tc- Parameters are /o = 10^ s“^ 
and KoM(O)^ = 1 eV. 


and vanishes at the Curie temperature Tc which also 
holds for AE{T) oc M{T)‘^. Thus, the system does not 
only switch faster with growing temperature due to an 
explicit increase of thermal fluctuations, but also due to 
a lower anisotropy barrier. This has a pronounced effect 
of the switching rate which is shown in Fig. 3. For conve¬ 
nience we use the decadic logarithm Ig = log^^g directly 
determine the exponent y of the frequency according to 
/ = 10^ s“^. Strong deviations from an Arrhenius law 
with a constant energy barrier are clearly visible. 

These deviations are not captured by a linear ad¬ 
justment of the Arrhenius law with a temperature shift 
T T— To, which yields to the Vogel-Fulcher law. Above 
Tc, this is obvious as the magnetization, and with it the 
energy barrier, vanishes. Close to Tc, the explicit values 
of /(T ~ Tc) loose their physical meaning since the con¬ 
cept of a switching rate with a separation of time scales 
breaks down [3]. However, from the existence of a finite 
critical Tc with AE{T > Tc) = 0 and f{T > Tc) = /o, 
we have to conclude that such a system can hardly be 
described by a constant energy barrier as in this case, 
/A£;=const(T) would approach /o for T ^ oo only. Thus, 
an accurate description of the switching rate requires the 
knowledge of the full temperature-dependent magnetiza¬ 
tion or even of the energy barrier at every temperature. 

The generalized Neel-Brown theory of the switching 
rate has profound consequences for the interpretation of 
experimental data. In particular, the interpretation of 
the prefactor in terms of an attempt frequency has to 
be revisited. In real systems, it is not a realistic task 
to measure the switching rate /(T) in the entire tem¬ 
perature regime from T = 0 to Tc- Usually, only a fi¬ 
nite (and small) temperature interval [Tmin,Tmax], with 
AT = Tmax — Tmin ~ 30... 50 K [7, 8] is experimen¬ 
tally accessible. Outside of this range, the switching is 
either too fast to be observable with present day tools or 


too slow to be measured in a reasonable time. In fact, 
AT is usually small enough such that the magnetization 
does not significantly change with temperature and the 
switching rate appears linearly depending on T~^ in this 
interval. Thus, it is tempting to fit the results to a linear 
function 

lg/ = lg/o®*-lg(e)^ (5) 

and to determine /o and Kq. However, the inter¬ 
pretation of the results for the fitted parameters /q^ 
and can actually be misleading. This is il¬ 

lustrated for a simple analytical example for which 
we assume for once a temperature-dependent magne¬ 
tization M(T) = M(0)y/l — TjTc. The switching 
rate resulting from Eq. (4) is now given by Ig/ = 
-lg/«* - lg(e)(ifoM(0)2)/(fcBT) with Ig/flt = Ig/o + 
lg(e)(iFoM(0)^)/(/cBTc'). This equation is in fact linearly 
depending on 1 /T and would actually yield the result for 
= iFoM(0)2. However, more importantly, the expo¬ 
nent of the fitted prefactor (which is often interpreted as 
the attempt frequency) differs from the physical one by 
lg(e)iFoM(0)^/(/i:BTc'). This difference can be very pro¬ 
nounced and, e.g., is for /o = 10^ s“^ and iFoM^(O) = 1 
eV given by Ig/o^ ^ Ig/o + . Hence, the attempt 

frequency /o = 10^ is in this case overestimated by five 
orders of magnitude for a typical Tc = lOOOK. 

For the three cases of M(T) used above, also the de¬ 
viation of the prefactor becomes now temperature de¬ 
pendent (this is not the case for the simplified model 
used in the previous paragraph). To show that this ef¬ 
fect can be very pronounced, we use the same model pa¬ 
rameters as above and numerically perform linear fits 
according to Eq. (5) in a small temperature interval 
T G [T^^ — Tc'/200,T^^]. The results are shown in Fig. 
4. The linear fits in general overestimate /o by several 
orders of magnitude. Interpreted as attempt frequencies, 
such linear fits yield to unrealistically large values re¬ 
ported in the analysis of experimental data on the basis 
of an Arrhenius or a Vogel-Fulcher law [6-9]. Naturally, 
the quality of the fit with respect to /o is the better, the 
less the magnetization depends on T in the respective 
window used for the fitting. Then, the magnetization 
is well approximated by the constant For a Bril- 

louin function, the (squared) magnetization varies less 
at lower temperatures (cf. Fig. 2). Thus, the best fits 
are achieved close to zero temperature. On the contrary, 
m(T)^ = [1 — (T/Tc)]^ changes the least close to the 
Curie temperature, which is also reflected in Fig. 4. Ob¬ 
viously, fits with the assumption of a constant energy bar¬ 
rier cannot deliver reliable results for iFo, when AE(T) 
is actually changing with T. Even though Kq probably 
is not misdetermined by orders of magnitude as it is fit¬ 
ted directly (and not its logarithm at T ^ oo as for /o), 
this correction of the barrier height has to be taken into 
account as well. 
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FIG. 4. (Color online) Linear fits (dashed lines) to the switch¬ 
ing rate Ig/ = Ig/o - lg(e)KoM(T)^/(A:BT) (full line), per¬ 
formed in a finite temperature interval around the given Tfit 
indicated by the vertical dotted lines (except for Tfit = O.lTc 
which lies out of bounds) for (a) a Brillouin function for S = 
4, (b) a Langevin function, and (c) a linear dependence (cf. 
Fig. 2). The exponents of /q^ follow for TcjT 0. Clearly, 
they can differ from the original /o by several orders of mag¬ 
nitude. The parameters are /o = 10^ s“^, Ko/M{0)^ = 1 eV 
and Tc = 500 K. 


The important consequences of a finite temperature 
window for the exponent of the fitted prefactor are il¬ 
lustrated in more detail in Fig. 5. We vary the posi¬ 
tion of the temperature window used for the fits and 
also the Curie temperature. In general, the regions for 
the optimal temperature are confirmed. An almost 
temperature-independent magnetization yields a well fit¬ 
ting attempt frequency /q^ ~ /o = 10^ Everywhere 
else, it is significantly overestimated. What is more, 
the lower the Curie temperature is, the worse are the 
fitted results, which can be again traced back to the 
temperature dependence of M(T). As the derivative 
dfM‘^{T) with respect to T = TjTc remains constant 
for all Tc for the three cases treated here, the derivative 
dTM‘^{T) = T^^dfM‘^{T) is proportional to the inverse 
Curie temperature. Thus, we find a larger deviation for 
a lower Tc. This is important, since superparamagnetic 
phenomena are typically investigated for very small sam¬ 
ples where the Curie temperature may be notably smaller 
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FIG. 5. (Color online) Fitted exponent of the attempt fre¬ 
quency. The magnetization is calculated by (a) a Brillouin 
function for S' = (b) the Langevin function (classical mag¬ 

netic moment) and (c) by a linear approach (cf. Fig. 2). The 
correct exponent would be lg/o(s“^) = 9. 


than in a bulk sample [18]. 

In this work we have shown that a temperature- 
dependent magnetic anisotropy in a superparamagnetic 
system can give rise to significant corrections to the 
standard Neel-Brown-Arrhenius law. For a correct in¬ 
terpretation of the parameters, the energy barrier for 
the switching of magnetic moment has to be based on 
a temperature-dependent magnetic anisotropy even in a 
simple model. This generates considerable corrections 
to the Arrhenius law. In particular the prefactor can no 
longer be determined in general by a fit of measured data 
on the basis of a temperature-independent anisotropy. 
Even if the fit appears plausible in a given finite temper¬ 
ature range, the attempt frequency can be misdetermined 
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by several orders of magnitude. 
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